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Background: Prostate cancer mortality has been decreasing in several high income countries and previous studies analysed the 
trends mostly according to geographical criteria. We aimed to identify patterns in the time trends of prostate cancer mortality 
across countries using a model-based approach. 

Methods: Model-based clustering was used to identify patterns of variation in prostate cancer mortality (1980-2010) across 37 European, 
five non-European high-income countries and four leading emerging economies. We characterised the patterns observed regarding the 
geographical distribution and gross national income of the countries, as well as the trends observed in mortality/incidence ratios. 

Results: We identified three clusters of countries with similar variation in prostate cancer mortality: pattern 1 ('no mortality 
decline'), characterised by a continued increase throughout the whole period; patterns 2 ('later mortality decline') and 3 ('earlier 
mortality decline') depict mortality declines, starting in the late and early 1990s, respectively. These clusters are also homogeneous 
regarding the variation in the prostate cancer mortality/incidence ratios, while are heterogeneous with reference to the 
geographical region of the countries and distribution of the gross national income. 

Conclusion: We provide a general model for the description and interpretation of the trends in prostate cancer mortality 
worldwide, based on three main patterns. 



Worldwide, prostate cancer is the second most frequent cancer 
and the sixth leading cause of oncological death among men. In 
2008, it accounted for approximately 250 thousand deaths (Ferlay 
et al, 2010) and 10.9 million of disability adjusted life years 
(Soerjomataram et al, 2012). 

Over the past two decades, the mortality rates have been 
declining in several high income countries (Bray et al, 2010; 
Bosetti et al, 2011). This has been explained in terms of improved 
early detection and treatment (Collin et al, 2008), although the 
effectiveness of prostate-specific antigen (PSA) screening (Andriole 
et al, 2012; Schroder et al, 2012) and radical prostatectomy remain 
open to discussion (Bill-Axelson et al, 2011; Wilt et al, 2012). 
Nevertheless, radiotherapy as primary treatment improves overall 
survival (Warde et al, 2011), immediate postoperative radiotherapy 
improves 5-year biochemical progression-free survival (Wiegel 
et al, 2009) and hormone therapy is associated to longer time to 



disease progression in metastatic disease and longer overall survival 
when associated with radiotherapy in locally advanced disease 
(Bolla et al, 2002). 

Previous studies grouped the countries mostly following 
geographical criteria (Hsing and Devesa, 2001; Oliver et al, 2001; 
Center et al, 2012) or according to trends after the introduction of 
widespread PSA testing (Bouchardy et al, 2008). Model-based 
clustering may allow a more meaningful grouping of the different 
settings with no a priori constraints, according to the mortality 
rates at onset of the observation period, as well as the magnitude 
and shape of its variation. 

We aimed to identify patterns in the time trends of prostate 
cancer mortality across countries using a model-based approach 
and characterised the patterns observed regarding the geographical 
distribution and gross national income (GNI) of the countries, as 
well as the trends observed in mortality/incidence ratios. 
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MATERIALS AND METHODS 



We analysed European countries, selected non-European high- 
income countries and selected leading emerging economies 
(United Nations, 2010) (Figure 1). 

Official prostate cancer mortality data were obtained from the 
World Health Organisation (WHO) database (World Health 
Organization, 2012) for the period between 1980, or the first 
calendar year with available data since then, and 2010, or the most 
recent year with data available, using the update of WHO database 
in 24 November 2011; countries with data available for < 10 



calendar years were excluded (Figure 1). Three different revisions 
of the International Classification of Diseases (ICD) were used to 
classify prostate cancer deaths in this period; we extracted the 
number of deaths corresponding to the codes A057 (ICD-8), B124 
(ICD-9) and C61 (ICD-10), as applicable. 

Estimates of the resident population in each mid-year were 
obtained from the 2010 revision of United Nations World 
Population Prospects (United Nations, 2010), even when mortality 
data refers to <90% of the country population (World Health 
Organization, 2012): Albania (70.8%), Brazil (79.7%), Republic of 
Moldova (83.2%), Serbia (83.1%), South Africa (76.9%) and The 
Former Yugoslav Republic of Macedonia (89.2%). 



MORTALITY DATA 

Source : WHO database (1980-2010) 
Age-standardised by world standard population 

European 

(ALB; AUT; BEL; BGR; BIH; BLR; CHE; CYP; CZE; DEU;DNK; ESP; 
EST; FIN; FRA; GBR; GRC; HRV; HUN; IRL; ISL; ITA; LTU; LUX; 
LVA; MDA; MKD; MLT; MNE; NLD; NOR; POL; PRT; ROM; RUS f ; 
SRB; SVK; SVN; SWE; UKR) 
Non-European more developed 
(AUS; CAN ;JPN; NZL; USA) 
Leading emerging economies 
(BRA; CHN*; IND; RUS f ; ZAF) 



INCIDENCE DATA 

Source : CI5plus database (1980-2002) 
Age-standardised by world standard population 

European 

(AUT; CHE; CZE; DEU; DNK; ESP; EST; FIN; 
FRA; GBR; ISL; ITA; NLD; NOR; POL; SVK; 
SVN; SWE) 
Non-European more developed 
(AUS; CAN; JPN; NZL; USA) 
Leading emerging economies 
(HKG) 



1 country excluded because there was no 
data available from the WHO database 
(IND) 



3 countries excluded because the data 
available in the WHO database referred to 
less than 10 calendar years 
(BIH; CYP; MNE) 





49 countries* 






f 


> 




48 countries 






f 


> 




45 countries 



25 countries 
(with incidence data collected) 



Pattern 1 
(18 countries) 
(ALB; BGR; BLR 
BRA; EST; GRC: 
HKG; HRV; JPN 
LVA; MDA; MKD: 
POL; ROM; RUS* 
SRB; SVN; UKR) 



Pattern 2 
(12 countries) 
(BEL; CZE; DNK; 
FIN; LTU; MLT; 
NOR; NZL; PRT; 
SVK; SWE; ZAF) 



Pattern 3 
(15 countries) 
(AUS; AUT; CAN; 
CHE; DEU; ESP; 
FRA; GBR; HUN; 
IRL; ISL; ITA; LUX; 
NLD; USA) 



Figure 1. Flowchart of the model-based approach used to identify prostate cancer mortality patterns and estimation of mortality/incidence 
ratios in the same countries. *For China, we used data from Hong Kong Special Administrative Region because Mainland China uses a special list 
of causes that does not encode prostate cancer. ''Russian Federation was considered both as an European country and a leading emerging 
economy; ALB = Albania; AUS = Australia; AUT = Austria; BEL= Belgium; BIH = Bosnia and Herzegovina; BGR = Bulgaria; BLR= Belarus; 
BRA= Brazil; CAN = Canada; CYP = Cyprus; CHE = Switzerland; CHN = China; CZE = Czech Republic; DEU = Germany; DNK=Denmark; 
ESP = Spain; EST = Estonia; FIN = Finland; FRA= France; GBR = United Kingdom; GRC = Greece; HKG = Hong Kong Special Administrative 
Region; HRV = Croatia; HUN = Hungary; IND = lndia; IRL=lreland; ISL=lceland; ITA= Italy; JPN = Japan; LTU = Lithuania; LUX = Luxembourg; 
LVA=Latvia; MDA= Republic of Moldova; MNE = Montenegro; MKD = The Former Yugoslav Republic of Macedonia; MLT=Malta; NLD = The 
Netherlands; NOR = Norway; NZL = New Zealand; POL = Poland; PRT = Portugal; ROM = Romania; RUS = Russian Federation; SRB = Serbia; 
SVK = Slovakia; SVN = Slovenia; SWE = Sweden; UKR=Ukraine; USA=United States of America; ZAF=South Africa. 
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Age-standardised mortality rates (ASMR), and the correspond- 
ing standard errors, were calculated by the direct method, using the 
world standard population (Doll and Smith, 1982) as reference, for 
the age groups, >45, 45-64 and ^65 years. 

Mixed models were used to describe the time trends in the 
ASMR. All models included random terms by country for the 
intercept, slope, quadratic and cubic terms. The observed mortality 
rates and the model predictions are presented in Appendix 1. 

Model-based clustering was used to identify groups of countries 
that share similar time trends in the ASMR, in the period 1980- 
2010, while distinguishing them from other homogeneous groups 
of countries regarding the variation in mortality rates. According 
to this method, the data (intercept, slope, quadratic and cubic 
terms) are assumed to have a multivariate normal distribution, 
parameterised by their means and covariances, generated by 
clusters. The geometric features (orientation, volume and shape) of 
the distributions are estimated from the data, and can be allowed to 
vary between clusters, or constrained to be the same for all clusters. 

In Appendix 2, we present the Bayesian Information Criterion 
(BIC) for models considering different geometric features of the 
distribution of the data and number of clusters. We considered the 
type of model and number of clusters corresponding to the lowest 
value of BIC as the most appropriate. 

The reliability of the model-based clustering was evaluated by 
10-fold cross validation (Efron and Tibshirani, 1997). The sample 
was divided into 10 partitions, and each of the subsets of 9 out of 
the 10 partitions was used to construct 10 different models. The 
agreement between the predictions from these models and those 
from the model constructed with the complete data set was 
calculated; the overall kappa coefficient was estimated by the mean 
of the kappa coefficients regarding the agreement between each of 
the 10 subset models and the full model. Data analysis was 
conducted using the software R 2.14.1. 

The patterns identified through the model-based approach were 
further characterised regarding the countries' region and GNI as 
well as the trends in mortality/incidence ratios. 

Age-standardised incidence rates (ASIR; world standard popu- 
lation) and the corresponding standard errors, for the age group 
^45 years, were abstracted from the Cancer Incidence in Five 
Continents database CI5plus for all the years with available data in 
the period 1980-2002 (International Agency for Research on 
Cancer, 2010). The incidence estimates reflect < 15% of the national 
population in nine countries (Austria, France, Italy, Japan, Nether- 
lands, Poland, Spain, Switzerland and United States of America 
(USA)) (Curado et al, 2007). Brazil and Germany were excluded 
from this analysis, because incidence data available refer to < 5% of 
the country population. We considered incidence data from 25 of 
the countries eligible for the analysis of mortality trends (Figure 1). 

To obtain a summary measure of the relation between the 
incidence and mortality rates, we computed ASMR/ASIR (MI) 
ratios for the countries and periods from which the necessary data 
were available (Appendices 3 and 4). For the United Kingdom (UK), 
incidence data were available only from Scotiand and England, and 
we computed the corresponding MI ratios using the specific 
incidence data and mortality data referring to the whole of UK. 



RESULTS 



We identified three clusters of countries with similar variation in 
prostate cancer mortality between 1980 and 2010, including 
pattern 1 (18 countries), pattern 2 (12 countries) and pattern 3 
(15 countries) (Figure 2 and Table 1). The reliability of the model- 
based clustering was good (kappa = 0.72). 

Pattern 1, labelled as 'no mortality decline', is characterised by 
the lowest rates in 1980, nearly half the observed for the other 
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Figure 2. Prostate cancer age-standardised mortality rates* (45 + 
years, direct method, World standard population), for each pattern 
identified. *Mean of the predictions for each of the countries included 
in the same pattern. 



patterns, and a steady increase between 1980 and 2010; the 
mortality rates almost doubled, from 24.2 out of 100 000 in 1980 to 
40.8 per 100 000 close to 2010. 

Patterns 2 and 3 are characterised by similar ASMR in 1980 
(close to 55 per 100 000 at age ^45 years), as well as by an 
inflection in the increasing trends between 1980 and 2010. In both 
patterns, the mortality rates declined approximately 40% from the 
highest values observed in this period to the rates in 2010. They 
differ essentially in the year when the decline starts; the highest 
rates observed in this period are observed nearly a decade earlier in 
pattern 3, and in the rates estimated for 2010, which are also lower 
for pattern 3. Patterns 2 and 3 were labelled 'later mortality decline' 
and 'earlier mortality decline' patterns, respectively. 

The main distinctive features of the patterns were essentially the 
same when only the age groups 45-64 or 65+ years were 
considered, although in the younger the inflection in the rates 
started sooner in all the patterns (Table 1). 

Although Southern and Eastern European countries are 
predominant in pattern 1 ('no mortality decline'), and Northern 
and Western European countries were assigned mainly to patterns 
2 ('later mortality decline') and 3 ('earlier mortality decline'), 
respectively, all clusters include countries from at least five of the 
six regions/groups considered (Figure 3). 

Regarding the distribution of the GNI per capita, pattern 1 
included mostly countries that were below the World Bank cutoff 
to define upper middle-income countries (3036 USD, in 1995 (The 
World Bank, 2012)) and pattern 3 included almost exclusively high 
income countries (>9385 USD, in 1995 (The World Bank, 2012)). 
However, the distribution of the GNI per capita in pattern 2 ('later 
mortality decline') overlapped the observed for patterns 1 and 3 
(Figure 4). 

In the 'earlier' and 'later' mortality decline patterns, the MI 
ratios were < 0.5 since the beginning and reached figures as low as 
0.2, which reflects the mortality declines, but mostly marked 
upward trends in the incidence rates; in many countries, the 
absolute changes in ASIR were higher than 100 per 100 000. 
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Table 1. Characterisation of the prostate cancer mortality patterns regarding the estimated age-standardised mortality rate (ASMR) (world standard 
population) in 1980, 1995 and 2010; percentage of changes in rates in the periods 1980-1995 and 1995-2010 and highest rates and corresponding year 
observed between 1980 and 2010, by age group 













Mortality rates (/1 00 000), median (P25-P75) 




Pattern 1 ('no mortality decline') 


Pattern 2 ('later mortality decline') 


Pattern 3 ('earlier mortality decline') 


45 + years 


ASMR a 


1980 
1995 
2010 


24.2 (17.4 to 31.4) 
31.6 (26.5 to 41.8) 
39.9 (36.1 to 50.7) 


54.7 (44.0 to 62.2) 
67.3 (57.1 to 72.0) 
47.6 (43.6 to 65.0) 


53.4 (50.5 to 57.3) 
62.2 (59.6 to 67.0) 
39.4 (34.5 to 46.7) 


Variation in ASMR (%) b 


1980-1995 
1995-2010 


38.3 (30.3 to 51.1) 
29.5 (15.2 to 50.0) 


20.3 (19.3 to 26.8) 
-22.6 (-36.5 to - 11.6) 


15.1 (8.2 to 19.0) 
-34.5 (-40.2 to -29.6) 


ASMR C 


Highest value 

Year of highest value 


40.8 (38.3 to 51.1) 
2007 (2002-2009) 


74.6 (66.5 to 81.3) 
1999 (1996-2001) 


68.1 (64.8 to 72.9) 
1992 (1991-1995) 


45-64 years 


ASMR a 








1980 
1995 
2010 


5.8 (5.2 to 7.1) 
8.6 (6.8 to 10.0) 

8.9 (7.0 to 11.7) 


9.4 (8.1 to 11.1) 
12.0 (10.7 to 14.0) 
8.7 (6.7 to 9.6) 


9.2 (8.4 to 9.8) 
10.3 (9.8 to 10.8) 

6.3 (5.3 to 7.6) 


Variation in ASMR (%) b 






1980-1995 
1995-2010 


43.8 (21.8 to 54.8) 
12.8 (-7.8 to 25.2) 


24.3 (15.1 to 41.0) 
-32.8 (-35.8 to - 17.1) 


10.8 (-1.1 to 21.1) 
-35.6 (-42.5 to -25.7) 


ASMR C 


Highest value 

Year of highest value 


11.6 (8.8 to 13.4) 
2002 (1998-2005) 


15.3 (13.2 to 17.0) 
1995 (1989-2000) 


13.2 (12.7 to 14.7) 
1990 (1988-1992) 


65 + years 


ASMR a 


1980 
1995 
2010 


72.3 (51.0 to 98.8) 

93.4 (76.5 to 129.4) 
125.7 (113.0 to 162.6) 


182.5 (144.2 to 200.4) 
219.9 (183.6 to 231.9) 
149.0 (140.3 to 218.1) 


178.5 (166.0 to 186.7) 
201.1 (196.0 to 218.6) 
127.3 (111.9 to 157.2) 


Variation in ASMR (%) b 


1980-1995 
1995-2010 


43.0 (29.4 to 50.7) 
34.7 (18.5 to 53.6) 


19.9 (18.8 to 30.7) 
-22.7 (-36.2 to -8.6) 


14.8 (8.8 to 17.9) 
-35.4 (-40.5 to -28.8) 


ASMR C 


Highest value 

Year of highest value 


126.3 (111.7 to 164.1) 
2007 (2002-2009) 


245.0 (214.2 to 255.2) 
2000 (1997-2000) 


224.4 (210.3 to 239.9) 
1993 (1991-1995) 


Abbreviations: ASMR = age-standardised mortality rate; P25 = percentile 25; P75 = percentile 75. 
a ASMR, estimated age-standardised mortality rate (model predictions). 

^Variation in ASMR (%), relative percentage of change in the estimated age-standardised mortality rate. 
C ASMR, observed age-standardised mortality rate. 



Regarding the trends in ASIR (Appendices 3 and 4), in most of 
the countries classified in pattern 3 ('earlier mortality decline') we 
can identify three periods with distinct patterns of variation in 
incidence and mortality. First, there is a steady increase in the 
ASIR, also along with increasing ASMR, though at a smaller pace 
than incidence, approximately up to the late 1980s or the early 
1990s. In a second period, usually lasting <5 years, the ASIR 



increased much more steeply than before while, the annual percent 
increase in ASMR remained similar to the observed in the first 
period. The third period is characterised by a decline in ASMR, 
while ASIR are still increasing. 

This is summarised in Figure 5, showing the trends in the 
MI ratios, which were above 0.4 in the early 1980s, decrease 
slightly in the first period, declined steeply in the second period, to 
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Pattern 1 Pattern 2 Pattern 3 



('no mortality decline')* ('later mortality decline') ('earlier mortality decline') 

■ Eastern Europe ■ Nortern Europe 

□ Southern Europe □ Western Europe 

a Non-European high-income countries » Leading emerging economies 

Figure 3. Number of countries from each region, by pattern of prostate cancer mortality trends. *ln pattern 1 , Russian Federation is represented 
twice, among Eastern European countries and the selected leading emerging economies. 



Pattern 1 
('no mortality decline')t 



Pattern 2 
('later mortality decline')}: 



Pattern 3 
('earlier mortality decline')¥ 




JPN 



HUN 




CHE 
LUX 



0 10 000 20 000 30 000 40 000 

Gross national income (USD per capita)* 

Figure 4. Gross national income per capita, by pattern of prostate cancer mortality trends. *The gross national income per capita (Atlas method) 
(GNI) in 1995 (the midpoint of the period under analysis) was obtained from the World Bank database (The World Bank, 2012). For Croatia, we 
used the GNI in 1997; for Serbia, we used the GNI in 1999; and for Estonia and Czech Republic, we used the GNI in 2002, as these were the first 
years with available data for these countries. According to the World Bank classification of economic development, this pattern includes two 
countries classified as low income (GNI <765 USD), eight countries classified as lower middle income (GNI between 766 and 3035 USD), five 
countries classified as upper middle income (GNI between 3036 and 9385 USD) and three countries classified as high income (GNI >9385 USD). 
'According to the World Bank classification of economic development, this pattern includes one country classified as lower middle income (GNI 
between 766 and 3035 USD), four countries classified as upper middle income (GNI between 3036 and 9385 USD) and seven countries classified 
as high income (GNI >9385 USD). "According to the World Bank classification of economic development, this pattern includes one country 
classified as upper middle income (GNI between 3036 and 9385 USD) and 14 countries classified as high income (GNI >9385 USD). 
CHE = Switzerland; GRC = Greece; HKG = Hong Kong Special Administrative Region; HUN = Hungary; JPN = Japan; LUX = Luxembourg; 
USD = United States dollars. 



<0.3, and still present a downward trend in the third period, 
to reach values close to 0.2 in 2002. Pattern 3 ('earlier mortality 
decline') presents the lowest values of MI ratio during the 
whole period of analysis, and it started to decline around 1985, 
approximately 7 years before the decline in mortality rates in this 
pattern. 

Pattern 2 (later mortality decline') is characterised by an almost 
stable MI ratio up to 1990, with values close to 0.5, and then by a 



decline until the end of the period, with values close to 0.30 in 
2002. The decrease of the MI ratio occurred about 9 years before 
the inflection of the upward trends in mortality rates, and 
approximately 5 years later than in the countries classified in 
pattern 3 ('earlier mortality decline'). Within pattern 2 ('later 
mortality decline'), Slovakia and Denmark are an exception to this 
general model of variation; in both the countries, the MI ratios 
were above 0.4 throughout the whole period and the ASMR 
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0.7 




0.1 



0.0 -I ! ! 1 1 

1980 1985 1990 1995 2000 

Year 

— Pattern 1 ('no mortality decline')t 
Pattern 2 ('later mortality decline')^ 

- - Pattern 3 ('earlier mortality decline')¥ 

Figure 5. Mortality/incidence ratio for prostate cancer* (age 45 + 
years), for each pattern identified, includes the countries for which it 
was possible to compute Ml ratios between 1988 and 2002 (Estonia, 
Hong Kong SAR, Japan, Latvia; Poland, Slovenia); 'includes the 
countries for which it was possible to compute Ml ratios between 1983 
and 2002 (Denmark; Finland; Lithuania; New Zealand; Norway; 
Sweden); 'includes the countries for which it was possible to compute 
Ml ratios between 1983 and 2002 (Australia; Canada; France; 
Netherlands; Iceland; Switzerland; United Kingdom*; United States of 
America). ^To represent the United Kingdom, we used incidence 
data from Scotland because incidence data available from England 
refers to less than the period considered for this analyses. *Each line 
represents the mean of the Ml ratios of the countries included in the 
same pattern. 

declined or stabilised, respectively, before steep increases in the 
ASIR could be observed. 

Pattern 1 ('no mortality decline') is characterised by a decrease 
in MI ratio since the early 1990s, although maintaining the highest 
values throughout the period of analysis. In 2002, the MI ratio was 
close to the observed in the countries classified in the other 
patterns when mortality started to decline. 



DISCUSSION 



Our analysis of recent trends in prostate cancer mortality identified 
three main patterns of variation in ASMR between 1980 and 2010. 
Pattern 1 is characterised by a continued increase in mortality 
throughout the whole period, while patterns 2 and 3 depict upward 
trends in the first years and subsequent declines, starting in the late 
and early 1990s, respectively. There was also a levelling of prostate 
cancer mortality worldwide, hence the difference between the main 
patterns was >50% in 1980 and <20% in 2010. 

We used an original approach to summarise the trends in 
prostate cancer mortality rates in different countries, and the 
unconstrained analysis of mortality patterns is the major strength 
of this paper. Model-based clustering has several statistical 
advantages over standard cluster approaches. It provides cluster 
solutions with 10 different covariance structures, while the most 
popular heuristic clustering methods (Ward's and K-means) only 
allow one type covariance structure (Fraley and Raftery, 2002), and 
also allows the choice of the model and the number of clusters to 
be recast as statistical model choice problems based on information 



criteria. The latter is also an advantage regarding factor analysis, 
where there is no model-based classification (Muthen, 2006). 

Our study adds to previous research on this topic the 
identification of three clusters of countries that are homogeneous 
regarding the variation in prostate cancer mortality, as well the 
trends in MI ratios, while remaining heterogeneous regarding the 
distribution of variables traditionally used for an a priori grouping 
of the countries, namely geographical region and GNI per capita. 

The MI ratio can be interpreted as a marker of the intensity of 
early diagnosis in different countries, as in the period under 
analysis its variation depended mostly on sharp increases in 
incidence, especially before mortality declines could be observed. 
Therefore, the assessment of its trends provides valuable empirical 
evidence for the discussion of the determinants of the patterns in 
mortality trends, given the scarcity of national data on the 
frequency of early diagnosis in several countries across time. 

In pattern 3 ('earlier mortality decline'), the steady increase in 
the ASIR even before the 'PSA era' probably reflects an increased 
diagnosis of prostate cancer due to the use of transurethral 
resection of the prostate to treat benign prostate hyperplasia (BPH) 
(Merril et al, 1999). The sharp increase around 1990 is consistent 
with widespread PSA testing. The estimates of the lead time 
associated with PSA screening range between 5 and 12 years 
(Draisma et al, 2003, 2009), but the lag between the steep increase 
in incidence rates and the mortality declines tends to be a few years 
only. 

Earlier prostate cancer detection coupled with the increased 
rates of radical prostatectomy for localised prostate cancer and the 
introduction of antiandrogen therapy that began before the 
introduction of PSA testing could contribute to the observed 
decline in ASMR (Feuer et al 1999; Mettlin, 2000). Moreover, 
during the 1980s and 1990s, there was a rapid rise in the use of 
radiotherapy and radical prostatectomy in the USA and in Europe 
for early disease (Selley et al, 1997), and in the 1990s there were 
improvements in those techniques (Peschel and Colberg, 2003). 

In pattern 2 ('later mortality decline'), the MI ratios were 
relatively stable up to 1990 and decreased steeply thereafter; at this 
point in time the MI ratios were also much higher in pattern 2 
(later mortality decline') than in pattern 3 ('earlier mortality 
decline'). The differences between the 'earlier' and later' mortality 
decline patterns may also be explained by the different extension to 
which early diagnosis is accomplished in each of these settings. In 
Denmark and Sweden, a steep increase in ASIR was seen 
approximately 5 years later than in the general trend of pattern 2 
(later mortality decline'). This could be the result of the limited use 
of PSA testing until around 1995 as a consequence of recommen- 
dations against opportunistic screening in asymptomatic men 
(Kvale et al, 2007), associated with the less frequent treatment for 
BPH (Moller, 2001). In Norway, the decline in ASIR after 2001 
corresponds to a decline in PSA testing, as a result of the issuing of 
recommendations to all general practitioners and urologists, 
similar to those made earlier in Denmark and Sweden 
(Norderhaug et al, 2003). 

Regarding pattern 1 ('no mortality decline'), Hong Kong Special 
Administrative Region (SAR), Japan, Latvia and Slovenia had a 
pronounced increase of incidence rates between 1990 and 1994, 
and in the remaining countries the ASIR increased at a steady rate. 
This pattern, however, is characterised by a gradual decline in the 
MI ratios, though never reaching values as low as those observed 
for the 'earlier' and later' mortality decline patterns. The inflection 
of mortality rates occurred only in Latvia and Japan, 12 and 15 
years after the start of the rapid increased of incidence, respectively, 
which is consistent with a late effect of PSA testing in these 
countries. 

The interpretation of our findings needs to account the accuracy 
of the mortality data. First, mortality rates were computed using 
the UN estimates for the overall population of each country and 
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the number of deaths in some countries refers only to a part of the 
population. However, the underestimation of the ASMR in those 
settings is not expected to compromise the comparison of the rates 
over time. Second, the adoption of three different revisions of ICD 
across the period considered could result in misclassification, 
although this is expected to have occurred in a small extent for 
prostate cancer (Klebba, 1980; Geran et al, 2005). Third, the 
reliability and validity of death certification may vary across 
settings and calendar years, especially in the selected leading 
emerging economies and some of middle-income countries 
included in our analyses. Although it is difficult to predict the 
extent to which these changes may have biased the magnitude and 
trends of the prostate cancer mortality rates, it is unlikely to have 
influenced the conclusions of our study. In fact, we identified three 
main patterns that are clearly distinct from each other and can 
accommodate differences between countries larger than could be 
expected from the heterogeneity in the quality of the data. This is 
supported by the fact that model-based clustering yielded 
essentially the same results when different subsets of the eligible 
countries were considered. However, some individual countries 
may have been misclassified regarding the cluster in which they 
were included, especially those that depict a pattern of variation 
compatible with more than one pattern. For example, South Africa 
was included in pattern 2 ('later mortality decline'), while 
presenting an ASMR in 1980 closer to the observed for the other 
countries classified in the same pattern, along with a variation of 
ASMR across time compatible with the observed for pattern 1 ('no 
mortality decline'). 

We opted for using incidence data obtained from the CI5plus 
database to ensure a higher homogeneity of data quality, which 
may be expected to minimise the potential bias associated with 
using data obtained from sources more heterogeneous regarding 
their quality standards. Therefore, the trends in the countries with 
sub-national data may not apply to the whole country, and the MI 
ratio estimates may be biased. However, our results and 
conclusions are robust and not dependent on accounting for any 
specific country; the trends in the MI ratios were essentially the 
same when sensitivity analyses (figures not shown) were conducted 
by excluding the countries with incidence data for <30% of the 
population or less than full coverage (though in the 'earlier 
mortality decline' pattern, only Iceland had incidence data for the 
whole population). 

Our analyses, including incidence data, could be also limited 
because the period covered by CI5plus extends only up to 2002. 
However, both for the 'earlier' and 'later' mortality decline 
patterns, no additional data is needed to interpret the mortality 
trends, as declines started to be observed before 2002, and the 
analysis of MI ratios in these countries sets a benchmark for the 
interpretation of the trends observed in the 'no mortality decline' 
pattern. Although no data were available for all the countries, the 
validity of our findings regarding the main objective of identifying 
patterns of variation in prostate cancer mortality is not 
compromised, as incidence data was used mostly for a finer 
characterisation and interpretation of the patterns identified. 

In conclusion, we provide a general model for the description 
and interpretation of the trends in prostate cancer mortality 
worldwide in the past decades, based on three main patterns 
homogeneous regarding the trends in mortality and MI ratio, each 
of them including countries from different regions and distinct 
categories of GNI. 
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Figure A1. Observed and predicted age-standardised mortality rate (direct method, world standard population), in the age group 45+ years. 
ALB = Albania; AUS = Australia; AUT = Austria; BEL= Belgium; BGR = Bulgaria; BLR= Belarus; BRA = Brazil; CAN = Canada; CHE = Switzerland; 
CHN = China; CZE = Czech Republic; DEU = Germany; DNK= Denmark; ESP = Spain; EST = Estonia; FIN = Finland; FRA=France; GBR=United 
Kingdom; GRC = Greece; HKG = Hong Kong Special Administrative Region; HRV = Croatia; HUN = Hungary; IRL=lreland; ISL=lceland; 
ITA= Italy; JPN = Japan; LTU = Lithuania; LUX = Luxembourg; LVA=Latvia; MDA=Republic of Moldova; MKD = The Former Yugoslav Republic 
of Macedonia; MLT= Malta; NLD = The Netherlands; NOR = Norway; NZL= New Zealand; POL= Poland; PRT= Portugal; ROM = Romania; 
RUS = Russian Federation; SRB = Serbia; SVK= Slovakia; SVN = Slovenia; SWE = Sweden; UKR = Ukraine; USA=United States of America; 
ZAF = South Africa. 
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Figure A2. Bayesian Information Criterion for models considering 
different geometric features (orientation, volume and shape) of the 
distribution of the data (intercept, slope, quadratic and cubic terms) 
and number of clusters. BIC — Bayesian Information Criterion (the plot 
depicts the BIC value multiplied by — 1);EII — Equal volume, equal 
shape, non-applicable orientation;VII — Variable volume, equal shape, 
non-applicable orientation;EEI — Equal volume, equal shape, 
coordinate axes orientation;VEI — Variable volume, equal shape, 
coordinate axes orientation;EVI — Equal volume, variable shape, 
coordinate axes orientation;WI — Variable volume, variable shape, 
coordinate axes orientation;EEE — Equal volume, equal shape, equal 
orientation;EEV — Equal volume, equal shape, variable 
orientation;VEV — Variable volume, equal shape, variable orientation; 
WV — Variable volume, variable shape, variable orientation. 
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APPENDIX 3 








Figure A3. Trends in age-standardised (45+ years, direct method, world standard population) mortality rates (ASMR), age-standardised (45 + 
years, direct method, world standard population) incidence rates (ASIR) and ASMR/ASIR (Ml) ratios. 
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Figure A3. Continued. 
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